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In this work we study the effects of the non-gravitational exchange energy ( Q ) between dark 
matter (DM) fluid and dark energy (DE) fluid on the background evolution of the cosmological 
parameters. A varying equation of state (EOS) parameter, oj, for DE is proposed. Considering an 
universe spatially flat, two distinct coupled models were examined to explore the main cosmological 
effects generated by the simultaneous reconstruction of Q and u> on the shape of the jerk parameter, j, 
through a slight enhancement or suppression of their amplitudes with respect to uncoupled scenarios, 
during its evolution from the past to the near future. In consequence, j could be used to distinguish 
any coupled DE models. Otherwise, the observational data were used to put stringent constraints on 
Q and oj, respectively. In such a way, we used our results as evidences to search possible deviations 
from the standard concordance model (ACDM), examining their predictions and improving our 
knowledge of the cosmic evolution of the universe. 
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I. INTRODUCTION 

The recent astronomical measurements of type la Su¬ 
pernovae Union 2.1 (Union 2.1 SNIa) composed by 580 
data M, the Baryon Acoustic Oscillation (BAO) de¬ 
tected in the clustering of the combined 2dF Galaxy Sur¬ 
vey (2dFGRS) and the Sloan Digital Sky Survey (SDSS) 
Data Release 7 (DR 7) main galaxy samples, the 6dF 
Galaxy Survey (6dFGS) and the WiggleZ Dark Energy 
Survey (WiggleZ) [H-Q, the observations of polariza¬ 
tion and anisotropies in the power spectrum of the Cos¬ 
mic Microwave Background (CMB: distance priors) data 
from the Wilkinson Microwave Anisotropy Probe 7 year 
(WMAP 7) fl0l - ITl , the observational Hubble (H) data 
set measured from galaxy surveys fl4l-fl8i]. and other, 
have confirmed that the present universe is undergoing an 
accelerated phase of expansion. In the literature, some 
theoretical approaches were taken into account to ex¬ 
plain this phenomenon, we are interested in an universe 
in where exists an exotic energy component with nega¬ 
tive pressure, named DE [19fi22l] , and which presumably 
began to dominate the evolution of the universe, only re¬ 
cently. Within this approach, the simplest candidate for 
DE is the Cosmological Constant A, which has an EOS 
parameter oj = —1.0 [23l - l26l |. Also, there exist other 
alternative models such as phantom model [271, quin- 
tom model [HJ, quintessence model )29[, the k-essence 
model Q, Chaplyging gas model [3l|, massive scalar 
field model Q and other. All these models predict dif¬ 
ferent dynamics of the universe. 

On the other hand, the properties of DE are mainly cha¬ 
racterized by oj. In such a way, due to our ignorance of 
its nature, it was parameterize empirically in a model in¬ 
dependent way. In this sense, we have followed two ways 
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to explore its behaviour. The first one was to paramete- 
rize oj in terms of some free parameters [33M371]. Among 
all the different parametrizations forms the Chevallier- 
Polarski-Linder (CPL) parametrization [34| is considered 
as the most popular ansatz oj = oj 0 + oj\[z/(\ + 2 )], 
where z is the redshift and oj 0 , oj\ are real parameters 
[34 ]. This ansatz has a divergence problem, when red- 
shift z approaches to —1 :38j. In addition, some non- 
parametric forms were found in [iiot . The second one 
was to choose an appropiated local basis representation 
for oj and after estimate the associated coefficients E3- 
[43l | . However, a divergence-free reconstruction for oj was 
proposed here, expanding oj in terms of the Chebyshev 
polynomials T n ,n £ N. To display how the method runs 
oj was expanded in terms of only the first three Cheby¬ 
shev polynomials T n , n = 0,1,2, and therefore, they are 
considered as a complete orthonormal basis on the finite 
interval [—1,1], and besides, belong to the Hilbert space 
L 2 of real values J44|. They were chosen because have 
the pro per ty to be the minimal approximately polynomi¬ 
als jl6l.l45|. 

On the other hand, within the universe another dark 
component has been assumed its existence, so-called 
DM, which acts exactly like the ordinary matter (pre¬ 
ssureless), but does not interact with DE, except gravi¬ 
tationally. The nature of these dark components are not 
still known and the possibility that within the universe 
exists a non-gravitational coupling in the dark sector 
could not be precluded, as well as, its possible effect on 
the dynamics evolution of the cosmological parameters 
should be considered 46 j59|. Some consequences of it 
were already studied in 6fl-63|, from which the strenght 
of the coupling should be very small. 

A huge amount of coupled DE models have already been 
investigated and fitted with cosmological data. Some 
of them were motivated by mathematical simplicity, for 
example, models in which Q oc Hp, in where H and 
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p denote the Hubble parameter and the energy density 
of dark sectors, respectively. It has three possibilities, 
namely, p = pdm (DM energy densit y), p = pde (DE 
energy density) and p = pdm + Pde [59| . On the con¬ 
trary, the models with Q oc p have been used in rehe¬ 
ating [if|, curvaton decay 0] and decay of DM into 
radiation [4S|. All these models strongly depend on the 
choice made for the Q form. So far, the coupled DE 
models have not been investigated, in a general form be¬ 
fore. In effect, some attempt to reconstruct Q from a 
general parametrization has been done by us in [64j] . In 
this paper we have considered different theoretical sce¬ 
narios, in where u> was always taken as constant. Here, 
we fixed Q oc H1q(1 + z) 3 , where the function Iq was 
reconstructed in terms of the first six Chebyshev poly¬ 
nomials. The analysis was done using a sample of SNIa 
Union 2 data. Our main results showed that the best 
fitted on Iq have prefered to cross the noncoupling line 
Iq = 0 during its evolution. 

The motivation of this article has been to go from the¬ 
ory to observations, following the prescription outlines 
by Cueva-Nucamendi 64(. Then, to follow this thread, 
a coupled model with two reconstructions describing to 
Q and w simultaneously has been proposed here. There¬ 
fore, we have postulated the existence of a general non- 
gravitational coupling between DE and DM (64] , intro¬ 
ducing a general phenomenological parametrization for Q 
into the equations of motion of these dark components. 
Here, Q was reconstructed expanding it in terms of the 
first three Chebyshev polynomials T n , n = 0,1,2. This 
has been the first attempt at reconstructing simultane¬ 
ously Q and ui from real data. 

Two distinct coupled DE models such as XCPL and DR 
were analysed here. Within these scenarios, the aim of 
our paper has been to study the effects that result from 
the reconstructions of Q and lo on the cosmological back¬ 
ground evolution, of some parameters such as (defined 
below) DM energy density parameter (Hum), decelera¬ 
tion parameter (q) and jerk parameter (j) [65H70|, whose 
amplitudes are modified with respect to those of uncou¬ 
pled models. Here, our models were constrained using an 
analysis combined of Union 2.1 SNIa [IH1], BAO jB-Q, 
CMB [Toi - ITB , and H data sets fl4l4l8j . 

Finally, we organize this paper as follows: The back¬ 
ground equation of motions for the energy densities, the 
definition of the geometrical parameters, and the recon¬ 
struction schemes for Q and ui are derived in section II. 
In section III we describe the coupled DE models worked. 
The priors considered and the observational constraints 
on the parameters space are discussed in section IV. We 
discuss our results in section V. In section VI we conclude 
our main results. 


II. BACKGROUND EQUATIONS OF MOTION 

In a flat Friedmann-Robertson-Walker (FRW) universe 
its background dynamics is described by the following set 


of equations for their energy densities (detailed calcula¬ 
tions are found in [6J], so we do not discuss them here.) 


Pb + 3 H pb — 0 , (1) 

p r + 4 H p r = 0 , (2) 

Pdm + 3 H p DM = +Q , (3) 

Pde + 3 (1 + w) H p DE = —<2 , (4) 


where pb, p r , Pdm and p EE are the energy densities of 
the baryon, radiation, DM and DE, respectively. Now 
defined the Hubble expansion rate as H = a/a, and also, 
indicates differentiation with respect to the time t. 

In what follows we shall assume that there is not ener¬ 
gy transfer from DE (DM) to baryon or radiation, and 
among them only exist a gravitational coupling 7l|. The 
critical densities p c = 3H 2 /8ttG, and the critical den¬ 
sity today p c , o = 3Hq/8ttG, in where Hq is the cur¬ 
rent value of the Hubble parameter, were conveniently 
defined. Considering that A = b,r , DM, DE, then the 
normalized densities are 


Pa _ Pa/ pc,o 


Ot = — = 

Pc Pc/ Pc,0 


o _ 

- -FTT > “40 = 


PA, 0 


E 2 Pc, o 

The first Friedmann equation is then given by 


E 2 = 


H 2 

m 


8t rG 


(Pb + Pr + PDM + pde ) , 

= [^6 + u* + n* DM + , 

and with the following relation for all time 
Qb + Q r + floM + ^de = 1 ■ 


( 5 ) 

( 6 ) 

( 7 ) 


The scale factor a is related with the redshift through 
a = 1/(1 + z), from which find dt/dz = —1/(1 + z)H(z). 
By substituting this last relation into Eqs. 0-0. and 
solving Eqs. 0-0, find the redshift evolution of 


n* b (z) 

n*(z) 

d PpM _ ^ ^DM 
dz 1 + z 
*DE _ 3(1 + U>) P k DE 
dz 1 + z 


fMl + z ) 3 , 
f^r,o(l + z) , 

~^DM^ Q(~) 

1 ~h z 

l + z 


( 8 ) 

( 9 ) 

( 10 ) 

( 11 ) 


Then, these equations have been fundamental to deter¬ 
mine the results within our models. 


A. Evolution of geometrical parameters 

The geometrical parameters of the universe are ob¬ 
tained by performing a Taylor series expansion of the 
scale factor a(t) around the current epoch, to- Conven¬ 
tionally, this series is truncated at a determined order 
[651] . Then, in this work we have been truncated such 
series at third order to study its behaviour, in where the 
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dimensionless coefficients such as deceleration parameter, 
q, and jerk parameter, j , are defined as (65| : 


q(z) s = - 1 + 


( 12 ) 


j( z ) = = Q( z ) + 2c i 2 ( z ) + (! + z )q'(z), (13) 

where d and a are the second and third derivatives of a 
with respect to time, respectively, and also, “ ' ” indicates 
differentiation with respect to z. Besides, some authors 
clehne the jerk parameter, j, with opposite sign [HE [63, 
Let us substitute Eqs. ©-Q into Eq. (THZ1) . 

q(z) = +-lo{z)£Ide{z) +--\ -(14) 

and its derivative is 

q'i z ) = +^u\z)n DE (z) + ^u(z)n , DE (z) + ^ : ^ . (15) 

These equations are frequently used in this work. 


III. DARK ENERGY MODELS 
A. ACDM model 

In this scenario, the function E 2 was found fixing both 
co(z) = — 1 and Q(z) = 0 into Eqs. (j5])- (fTTl) 


E 2 — ^b(z) + 11 *(z ) + Hdm,o (1 + z) 6 + £Ide, o 


( 20 ) 


moreover, q and q’ are given by the Eqs. m-m- then 


«(-') = +j 


b(z) + Qdm(z) + 2 fi r (z) — 1 , 

Q' b (z) + Q'dm( z ) +2 Q' r (z) , (21) 

from these equations, reconstructed the parameter j(z). 


= + 2 


B. CPL model 

B. Parametrizations of Q and w 


The Chebysliev polynomials form a complete set of or¬ 
thonormal functions on the interval [—1,1] and have the 
property to be the minimal approximating polynomials, 
which means that have the smallest maximum deviation 
from the true function at any given order EE [si. 

In general, an energy exchange is described as the cou¬ 
pling between both dak fluids and phenomenologically it 
is chosen as a rate proportional to H 

Q = Hp DM lQ ■ (16) 

Here, the strength of the coupling is characterized by Iq, 

2 

= ^Tu , (17) 

n =0 

in where the coefficients of the pol yno mial expansion A„ 
are free dimensionless parameters [64J and 

T 0 (z) = l, T\(z) = z , T 2 (z) = (2z 2 — 1) , (18) 

represent the first three Chebyshev polynomials. 

Within the CPL model, the past evolution history may be 
successfully described by its EOS parameter, w, but the 
future evolution may not be explained, because to grows 
increasingly, and then, encounters a divergence when z —> 
— 1. This is not a physical feature. Therefore, a novel 
reconstruction form for ui has been proposed here to avoid 
such divergence problem. Hence, to is defined as 

2 

U = ^ UJmTm , (19) 

771=0 

where w 0 ,wi and lo 2 are free dimensionless parameters. 
The Chebyshev polynomials of order m = 2 were defined 
by Eq. m- Using numerical simulations we will com¬ 
pute the best fitted values for Ao, Ai, A 2 , w Q , lo\ and ui 2 , 
respectively. 


Within this model, E 2 was determined replacing both 
w(z) = w 0 + vui [ 2/(1 + 2 )], where w 0 , w\ are real param¬ 
eters and Q(z) = 0 into Eqs. ([5])- (HTT) 


E 2 = 


S^6,o(l + z ) 3 + H rj o(l + 2 ) + Qdm, o(l + 2 ) 3 


+n DEfi (l + zf 1+w ° + "’ 1 >exp( 1 ^) 


• ( 22 ) 


Then, the following relation 


u\z) 


Wl 

(I + 2) 2 ’ 


(23) 


is substituted into Eqs. m-m and m, from which 
reconstructed q(z) and j(z), respectively. 


C. XCPL model 


Here, firstly a coupled model was defined putting both 
ui = uj 0 + wi(*/l + 2 ), where w D , are real free parame¬ 
ters and Q{z) given by Eqs. (HBl) - (fI71) . into Eqs. (l81)- ([TTjl . 
The explicit form for £l* DM and Sl* DE are reached, solving 
Eqs. m-m, respectively, 


dm ( 2 ) — (1 + z) 3 Q E M, osxp 


Zmax 

~2 


2 1 
^ ^ ^nln{ z ) > 

71=0 


VdeM = (1 + *) 3(1+,3i) 


^DB, 0 exp 



Zmax 

2 


^dm, oexp 


/ 3wi 

Vl + 2 


2 

T. A n S n (z,w) 

71=0 


(24) 


(25) 
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The following average integrals have been defined 

r r n (x) 

Jo (! + z) 


(1 + J)( 1 +3^) 


in where we have also defined the following expressions 
for all n £ [0,2] (see Appendix „ and HI) 


A(x) = exp y T Q3: y^A„I„(x)j 
A(x) = exp ( ~™ ax 

V n=0 / 

, {—3wi\ vx, , ( —3wi 

B(x) = exp ( --- ) , B(x) = exp 


(28) 


1+5 

r ( \ _ f X T n (x)'~ 

In \ z ) = / ?— rr~\ dx > 

J -1 {a + bx) 

O , r T n (x)A(x)B(x) ^ 

S " (Z ’"' ) = 1, (a + te)»+»> dx ' 


a + bx 


where the following relations were defined 

2 


dx « Zmax I (A 

UJy ■ — 2 i nw ) 1 

(26) 

C(z) 

1 ~ Zmax ri / — \ 

dx « o«(z, w) , 

(27) 

D{z) 



F(x) 


C(z) = SIde, o(l + ^) 3 exp 


3 Zr> 




n —0 


-s(hi (i + 5=[i + x]) 


Z max . .A' 


G{x) = 


z max , 

1 H-—(1 + X) 


1 —3(c*Jl—W2 —UJo) 


L(z) = exp 


3 (u>2Z + (wi — 2 w 2)-2 + Wi + W 2 ) 


Within this model, the function E 2 was constructed from 
Eqs. ©-©, ([H, (Hi) and ©. Furthermore, the rela¬ 
tion 


uj'(z) = Wi + 4ztc>2 , 


(30) 


is replaced into Eqs. dHD-dUD and m to reconstruct 
q{z) and j{z), respectively. The basic analytical expre¬ 
ssions for I n (x) and I n {x ) are found in Appendix D 


and the quantities, 


x — 2(z/z max ) 1 , 

, Zmax 

a = 1 -I-— , 


OJ = L0 o + +L , 


6 = 


Zmax 

2 


where z max is the maximum value of z in which the obser¬ 
vations are possible so that x £ [—1,1] and \T n (x)\ < 1. 
Therefore, the function E 2 was constructed from Eqs. 
©-©, (El) - G5l) and ©. Similarly, Eq. El) is then 
substituted into Eqs. tin-cog) and tin to reconstruct 
q{z) and j(z), respectively. 


D. DR model 


IV. CURRENT OBSERVATIONAL DATA AND 
COSMOLOGICAL CONSTRAINTS. 

In this section, we describe how we use the cosmo¬ 
logical data currently available to test and constrain the 
parameter space of our models proposed. 


A. Type la supernovae data set. 

For the SNIa observations, we consider “The Super¬ 
nova Cosmology Project” Union 2.1 composed of 580 
SNIa data. The distance modulus /i(z,X), is defined 
as the difference between the apparent ( m ) and absolute 
(M) magnitudes, so that their observed and theoretical 
values are 


Secondly a coupled model was modeled setting both 
w(z) = w 0 + wiz + it> 2 ( 2 z 2 — 1), where w 0 , to 1 , 0 J 2 are 
real parameters and Q(z) given by Eqs. (fTHl) - (fI71) . into 
Eqs. m -m- The explicit form for and £E de 

are reached, solving Eqs. tiua-«nj, respectively. In this 
way, Q and w were simultaneously reconstructed. For 
this model Eq. (El) represents the solution of Eq. 113, 
and hence, the solution of Eq. m was obtained using 
Eq. (l28l) and Appendix,, 


/i obs = m obs - M 


M th (z, X) = 51og 10 


Dl{z, X) 
Mpc 


+ MO ! 


(31) 

(32) 


where /io = 25 — logjn-ffo, and the Hubble-free the lumi¬ 
nosity distance D E [l[ in a flat cosmology is 


^.X) = (l+-sf (33) 


^*de( z ) = c {z)+D{z)L{z) y^uj n T n (x) 


n —0 


A(x)F(x) 

G(x) 


dx 


(29) 


in where Ff(z,X) is the Hubble parameter, i.e., Eq. © 
and, in general, X represents the model parameters 

X = (Hq, £U,0j Udm.o+o! wi, W 2 , Aq, Ai, A 2 ). (34) 
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The best fitting values of the parameters in a model are 
determined by the likelihood analysis of 


2 [/^ a (z fc )-/X th (s fc ,X)]‘ 


Xsn(X, to o) = 


fc=1 




(35) 


in where cr(zk) is the corresponding lcr error of distance 
modulus for each supernovae. The parameter po is a 
nuisance parameter. According to [|J, y 2 is expanded as 

Xs N (X) = A,- 2fj>oBi + Mo 2 Ci , (36) 


where 

1 h ^ a (^) 

580 ^o bs{zk) _ ^th {zkjX) } 

1 

580 

Cl 2 ' 


The Eq. (136l) has a minimum for po = B\jC\ at 


(37) 


xl N (X) = A 1 (X)-®| . (38) 

Since x|iv,mm = xliv,mm> instead minimizing x|jv we 
will minimize Xsatj which is independent of po- 


B. BAO data sets 


Eisenstein et al.@, first found a well-detected peak 
of the imprint of the recombination-epoch acoustic 
oscillations in the large-scale correlation function at 
100ft -1 Mpc (ft = Ho/100Kms~ 1 Mpc~ 1 ) separation 
measured from a spectroscopic sample of 46, 748 lumi¬ 
nous red galaxies of the SDSS. Also, Percival et al.Q, 
investigated the clustering of galaxies within the spec¬ 
troscopic SDSS-DR7 galaxy sample including, the lumi¬ 
nous red galaxy, main samples, and also the 2-degree 
Field Galaxy Redshift Survey (2dFGRS) data (in total 
893, 319 galaxies) observed BAO in power spectrum of 
matter fluctuations after the epoch of recombination on 
large scales. This allowed to detect the BAO signal at 
2 = 0.2 and z = 0.35. 

Eisenstein first and Percival after constructed an effective 
distance ratio D v (z), which encodes the visual distortion 
of a spherical object due to the non-Euclidianity of a 
FRW spacetime, defined as 


D v {z,X) = -±- 
B o 

[ (1 + *W(,)^)] 

1/3 

5 

c 

\( r dz> V 

Z 

Ho 

\Jo E(z',X)J E(z, X) 


1/3 


(39) 


z 

d z 

0.20 

0.1905 + 0.0061 

0.35 

0.1097 + 0.0036 


TABLE I. Summary of SDSS-2dFGS BAO data set 0. 


where D A (z) is the proper (not comoving) angular diame¬ 
ter distance, which has the following definition 


D a {z, X) = 


dz' 


(1 + z)J o H(z', X) ■ 

The comoving sound horizon size is defined by 

3 (a')da' 


,s_ rcs( 
M)=C\ —2 

Jo a' 


(40) 


(41) 


jo « 2 H{a') ’ 
being c s {a) the sound speed of the photon-baryon fluid 

1 


r 2 ( )= 6 _Z = }_ 

Cs[a> 5p 3 [1 + (30b/4f2 r )a 


Considering the Eqs. sid and (H21) for a z given 

rVCH-*)_ da 

0 


(42) 


r ' w = 7S, 


a 2 if(a)v/l + (3ft M /4fi 70 )a ’ 

(43) 

where flb,o and fl 7 ,o are the present-day baryon and pho¬ 
ton density parameters, respectively. In this paper, we 
have fixed fl 7j o = 2.469 x 10 -5 ft -2 , = 0.02246ft -2 , 

given by the WMAP7 data 0 and O r! o = 0 7j o(l + 
0.2271N e ff) with N e ff the effective number of neutrino 
species (here the standard value, N e ff = 3.04 was chosen 
fl3j|h The peak position of the BAO depends on the ra¬ 
tio of D v {z ) to the sound horizont size at the drag epoch 
(where baryons were released from photons) z d , which 
was obtained by using a fitting formula @]: 


1291(fl M ,0* 2 ) 0 ' 251 /i i l Jr~i u 

Zd — -i , r\ 1 9\n HTS (1 4“ ^1 (^6,0^ 


2\b 2 ) 


1 + 0.659(fl M ,o^ 2 ) 0 - 828 

where Qm,o = 9dm,o + I2&,o and 

fti = 0.313(f2M,o^ 2 ) -0 ' 419 [1 + 0.607(fl M ,o/i 2 )°- 674 ' 
b 2 = 0.238(fi M ,oft 2 )" 


. (44) 


2\0.223 


The distance radio d z at z = 0.2 and z = 0.35 extracted 
from the SDSS and 2dFGRS @-0, are listed in Table [1] 


do.2o(X) — 


r 8 (zd) 


D v { 0.2,X) ’ 


do. 35 (X) = 


r s (z d ) 


D v ( 0.35, X) ’ 

. (45) 

where r s (zj,X) is the comoving sound horizont size at 
the baryon drag epoch. Using the data of Table 0 then 
the inverse covariance matrix of BAO 0 is 


C -1 

^SDSS 


+30124 — 17227 \ 
-17227 + 86977 ) ' 


(46) 
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Z eff 

A(z) 

0.44 

0.474 ± 0.034 

0.60 

0.442 ± 0.020 

0.73 

0.424 + 0.021 


TABLE II. Summary of the effective redshifts z e // and the 
corresponding measurements for the acoustic parameter A(z) 
of the WiggleZ BAO data [i|. 


The errors of the SDSS-2dFGS BAO data are contained 
in CgQgg. The y 2 °f the SDSS-2dFGS BAO data set is 


xIdss(X) — (Adi) (CsDss)ij , 

where Adi = d‘^(X) — d° bs is a column vector 


df(X)-df s 


d 0 . 2 o(X) - 0.1905 \ 
do. 35 (X)- 0.1097 ) ’ 


(47) 


(48) 


and “t” denotes its transpose. 

On the other hand, in the low redshift region (z < 1), 
Beutler et al. studied the large-scale correlation func¬ 
tion of the 6dFGS and detected a BAO signal [§]. This 
detection allowed to constrain the distance-redshift re¬ 
lation D v (z ) at z e ff = 0.106, D v (z e ff) = 456 ± 27 
Mpc and a measurement of the distance ratio, d Zeff = 
r s {zd)/D v {z e ff ) = 0.336 ± 0.015, where r s (zd ) is the 
sound horizon at the drag epoch Zd- 
The x 2 function of the 6dFGS BAO data set is given by 


X6dFGSl 


/ d z - 0.336\ 
V 0.015 ) 


(49) 


Also, in Q, Blake et al. presented new measurements of 
the BAO at z = 0.44,0.6, and 0.73 in the galaxy corre¬ 
lation function of the final data set of the WiggleZ. This 
sample included a total of N = 158, 741 galaxies in the 
range 0.2 < z < 1.0. Here, the acoustic parameter A, 
was defined as 


A{z) = D v (z)y/Q^(^j , 



Using the BAO data showed into Table |TT1 the inverse 
covariance matrix of WiggleZ BAO Q is given by 


Ia{z») 

302.09 ± 0.76 

R{z t ) 

1.725 + 0.018 

z* 

1091.3 + 0.91 


TABLE III. The distance priors given by fl3j]. 


and “t” denotes its transpose. 

Since Eqs. C71) . fl^l) and O l. the total y 2 f° r the BAO 
data sets is constructed 

XBAO = XSDSS + XfidFGS + X WiggleZ ■ (54) 


C. CMB data set 


The Union 2.1 SNIa and BAO data sets contain in¬ 
formation about the universe at low redshifts, we now 
include WMAP 7 data [l3| to probe the entire expan¬ 
sion history up to the last scattering surface. The shift 
parameter R is provided by [11] 


Hn 


R(z*, X) = —a/Um,o(1 + z*)Da{z*,X.), 
c 


= ■s/A 


M, 0 


: dy 

E (y ) 


(55) 


where the distance Da and E(y) are given by Eqs. (TTO1) 
and m, respectively. Then, the redshift z » (the decou¬ 
pling epoch of photons) was obtained using the following 
fitting function 0 


z* = 1048 


1+0.00124(Q 6i0 /U) 


2\ —0.738 


\+gi(p,M,oh 2 ) 92 


(56) 


where 9 m,o = 9jjm,o + ^&,o, moreover, gi and g 2 are 


_ 0.0783 (n bfi h 2 )-°- 238 0.560 

91 ~ 1 + 39.5(U 60 h 2 ) 0 - 763 ’ 92 ~ 1 + 21.1 (Clboh 2 ) 1 - 81 ' 

(57) 

An angular scale l a for the sound horizon at decoupling 
epoch was defined as 


U(X) = (1 + z*) 


7 tD a (z*,X) 
r s {z*, X) 


(58) 


/ +1040.3- 807.5 + 336.8 \ 

^WiggleZ = -807.5 + 3720.3- 1551.9 . (51) 

\ +336.8 - 1551.9 + 2914.9 ) 

Then, define the y 2 function of the WiggleZ BAO data 

XwiggieZ(X) = (AAj) (c w ' iggleZ ) (A Ajf , (52) 

where A A* = Af 1 (X) — A° hs is a column vector 
/ A(0.44,X) - 0.474 \ 

Af(X)-A° bs = A(0.60,X)- 0.442 , (53) 

\ A(0.73,X) - 0.424 ) 


where r s (z*, X) is the comoving sound horizon at z*, and 
is given by Eq. (14311 . The maximum likelihood values 
according WMAP 7 data [l3j are given in Tabic IIIII 
Then, following [l3| the y 2 for the CMB data is 

Xcmb(X) = (Aa*) (C c ^ B ) y . (A xff , (59) 

where A Xi = xf h (X) — x° bs is a column vector 

/ l A (z*) ~ 302.09 \ 

xf l (X) - xf s = R(z*) — 1.725 , 

\ z* - 1091.3 J 


(60) 
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“t” denotes its transpose and (C c ^ B )y is the inverse 
covariance matrix ([131]) given by 


fi —i = 

l 'CMB — 


/ +2.3050 
+29.698 
\ —1.3330 


+29.6980 

+6825.27 

-113.180 


-1.3330 \ 
-113.18 
+3.4140 ) 


(61) 


The errors for the CMB data are contained in C c ^ B . 


D. Observational Hubble data (H) 


in UMi, the authors established that it is possible 
to compute the observational Hubble data, by using di¬ 
fferential ages of galaxies through the measuring of dz/dt. 
Then, the Hubble parameter was expressed in terms of 
the differential ages as 


H(*) = 


1 dz 
1 + z dt 


(62) 


Simon et al.[l6j found Hubble data over the redshift 
range [0,1.8]. In 0 , the authors found new data of 
the Hubble parameter at z £ [0.35,1] from SPICES 
and VVDS galaxy surveys, respectively, which are listed 
in Table M\ In addition, in 18j, the authors took the 
BAO scale as a standard ruler in the radial direction 
and found three more additional data: H(z = 0.24) 
= 79.69 ± 2.32, H(z = 0.34) = 83.80 ± 2.96 and H{z = 
0.43) = 86.45 ± 3.27 (in units of Kms~ 1 AIpc~ 

The x 2 f° r the observational Hubble data is [72| 


15 


Xh 


<*) = £ 


[g‘ t (X,Z„)-g° 1 '*(2 i )] ; 

<J 2 {Zi) 


(63) 


where X represents the parameters of the model, H th is 
the theoretical value for the Hubble parameter, H obs is 
the observed value, and er(+) is the standard deviation 
measurement uncertainty. Here the summation is over 
the 15 observational Hubble data at Zj. This test has 
been already used to constrain several models in [73l - l75| . 
Therefore, to fit our models with observations, we use all 
the data sets described above. 

The best fitted parameters are obtained by minimizing 

X — XsN + X BAO + x CMB + X H ■ (64) 

From (l64l) , we will construct the total probability density 
function pdf as 


pdf(X) = Ae * / 2 . (65) 

where A is a integration constant. 


E. Constant Priors 

In this work, we have assumed that baryonic matter 
( b) and radiation (r) are not coupled to DE or DM and 


are separately conserved 0. In this regard, we believe 
that the intensity of the interaction, Iq, is not affected 
by the values of and f2 ri o, respectively. Due to it, in 
this paper, we have fixed: f2 7i o = 2.469 x 10 _5 /i” 2 and 
flfj.o = 0.02246 h~ 2 , given by WMAP 7 data [Tj. Then, 
based in these assumptions, construct pdf for each of our 
models. The priors on the parameters space are given in 
Table El which were used in all the observational tests of 
our models. They have allowed us to compute the best 
fitting values of the free parameters. 


V. RESULTS 

In this section, we present the results of the fitting 
on the models listed in Table ED using the Union 2.1 
SNIa data set, the BAO data set, the CMB data from 
WMAP 7, the H data set and the priors described in Ta- 
bleEl Likewise, for the uncoupled ACDM and CPL mod¬ 
els, the corresponding free parameters to be estimated 
are: X = (SIdm,o,H 0 ) and X = Ho, w Q , a+). 

Meanwhile, for the coupled XCPL and DR models the 
free parameters are: X = (Ddm,o, Ho, u 0 , uj\, Ao, Ai, A 2 ), 
X= (n DM ,o,H 0 ,u} 0 ,ui,uj 2 , A 0 , Ai, A 2 ), respectively. 

In each model, the function Xr nm ■ the one-dimension 
probability contours, the best fitting parameters, and 
their errors at la and 2a were computed, by using the 
Bayesian statistic method, as are shown in Figures [l] [2] 
and [3] respectively. 

The values of the functions q , j. Iq, w and CIdm evalua¬ 
ted in z = 0 (today) are denoted as q 0 , j 0 , I Q , uj 0 and 
f2DM.o, respectively, and are presented in Table 1X171 
Furthermore, the coupled models show that the uj = — 1 
crossing feature is more favored by the reconstructed u> in 
the DR model with two crossings in the past (z = 21.11 
and z = 7895) than that found in the XCPL model with 
only one crossing in a recent epoch (z = 0.07) [fb • These 
crossing points obtained from the best reconstructed to 
are illustrated in Figure [D 

Let us now see Figure [D within the coupled models have 
considered that I_|_ denotes an energy transfer from DE 
to DM , instead, I_ denotes an energy transfer from DM 
to DE. A change of sign on the best reconstructed Iq is 
linked to the crossing of the noncoupling line Iq(z) = 0. 
In this regard, within the coupled models have found a 
change from I + in the past to I_ in the present and vice 
versa. According to this Figure and Table ETD note that 
a non-negligible value of Io at ler error has been found in 
the coupled models, and whose order of magnitude is in 
agreement with the results obtained in [55li59l 1641] . 

Due to the two minimums obtained in each coupled 
model (see Table ED, then two different possibilities to 
reconstruct Iq have been found here. Therefore, for any 
z range, an interesting mixture of I_|_ and I_ may be des¬ 
cribed there. 

In what follows we compared both the results of the 
XCPL model with those of the CPL model, and also, 
the predictions of the DR model with the corresponding 




z 

H(z) 

la 

0.0 

74.2 

± 3.6 

0.10 

69 

±12 

0.17 

83 

±8 

0.27 

77 

±14 

0.40 

95 

±17 

0.48 

97 

±60 

0.88 

90 

±40 

0.90 

117 

±23 

1.30 

168 

±17 

1.43 

177 

±18 

1.53 

140 

±14 

1.73 

202 

±40 


Parameters 

Constant Priors 

Qdm, 0 

Ho^kms^ 1 Mpc -1 ) 
Ao 

Ai 

A 2 

UJo 

Ui 

UJ2 

[0, 0.7] 

[20,140] 

[-1.5 x 10 +6 , +1.5 x 10+ 6 ] 
[-1.5 x 10 +6 , +1.5 x 10+ 6 ] 
[-1.5 x 10 +6 , +1.5 x 10+ 6 ] 
[-2.0, -0.3] 

[-1.0,+1.0] 

[-10,+10] 


TABLE V. Shows the priors on the parameter space. 


TABLE IV. Shows the observational H(z) data T7j] . 



fl DM,0 






fl DM.O 


FIG. 1. (color online) Shows the one-dimension probability contours for each of the parameters of our noncoupled models. 
In each panel the black star denotes the best fitting value of the parameters. The ler and 2er represent the errors. 


of the ACDM model. Otherwise, according to the results 
presented in Figures [I] and [5] For 2 > 3.98 the values of 
the amplitudes of Q dm in the coupled models are slightly 
modified by the values of Iq (1+ or I_) when they are 
compared with the uncoupled models. For I + > 0, the 
amplitudes of f l dm are amplified, instead, for I_ < 0, 
these amplitudes are sup pressed. These results coincide 
with those found in 1531 ] . In addition, in any another z 
region, these effects are not possible. They are described 
in Table rVTffi (left and center above table), in where the 
2 ranges and the coupled models are indicated. Likewise, 
we have also noted that, the shape of £Ide and the valu¬ 
es of its amplitudes are not significantly affected by the 
reconstructions of Iq and ui with respect to uncoupled 
models. Furthermore, we also confirm that the coinci¬ 
dence problem is alleviated in these coupled models, but 
they may not solve it. The below panels in Figure [5] as 
well as, the right above panel and below panels in Figures 
[G] and [TJ the constraints at la and 2er on VLdm , Udei Iq, 
to, q and j have been omitted to obtain a better visuali¬ 


zation of these effects. 

From Figure [5] (left above panel) our coupled models have 
predicted that a transition from a deceleration era at 
early times to an acceleration era at late times has been 
reached by the universe and the redshift for this change, 
Zt, is found to be zt ~ 0.75. 

Then, at z —> — 1 the DR model cannot determine, if 
the big-rip [77j may or may not occur in the universe, 
instead, for the XCPL model the universe will finish in 
a big-rip. These results coincide with those obtained in 
[ 3 ! ]. However, from Figure |5] (left above panel) and Table 
Em] (right above table) our results in the XCPL model 
have revealed that in determined z intervals the ampli¬ 
tude of q is slightly enhanced due to the increasing of the 
magnitudes of Iq and ui (see right above panel and be¬ 
low panels in Figure [6]) with respect to that of the CPL 
model. Such a effect was realized when Qdm becomes 
less concentrated at log(z) £ [—2, +0.3] (see left panel in 
Figure [5]). It was the epoch of the DE dominance, from 
which, the universe was led to an accelerated expansion. 








































9 






- 0.2 - 0.1 0 0.1 0.2 



FIG. 2. (color online) Displays the one-dimension probability contours for the parameters of the XCPL model and their 
constraints at la and 2cr, respectively. In each panel the black star denotes the best fitting parameter. 
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FIG. 3. (color online) Displays the one-dimension probability contours for the parameters of the DR model and their constraints 
at lcr and 2cr, respectively. In each panel the black star denotes the best fitting parameter. 


Regarding Figure G] (left above panel), we have found 
qualitatively different asymptotic values in the near fu¬ 
ture (z —> —1), but similar asymptotic values in the past 
(z 1) are exhibited by the best reconstructed j. Now 
let us analyze this Figure (right upper panel and lower 
panels) and Table fVTTTl (below tables), from which, note 
that the magnitudes of Iq and ui have imprinted new 
physical effects on the amplitudes of the parameter j. 
Within the coupled models the amplitudes of j were pro¬ 
gressively increased or reduced in determined z regions, 
with respect to those of the uncoupled models for in¬ 
creasing of the magnitudes of Iq and ui. Furthermore, 


the expansion of the universe was modified by the con¬ 
sequent diminution of Qdm- Therefore, we have shown 
that las magnitudes of Iq and u> are strongly related with 
the magnitudes of ft dm, q and j, respectively, as are seen 
in F igures [5j [6] and [7J 

We now compare our results with those obtained by other 
researchers. In [bit several parametrizations for ui were 
proposed, and also, a series for uj in [6(| (see Tables 1 
and 2), was investigated, and then, compared it with a 
expansion of the scale factor up the fifth order. The re¬ 
sults of [65] and @ are comparable with our results at 
lu error (see Tabic IVTIl) . Furthemore, in [67] and [6.8] 
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TABLE VI. Shows the best fitting cosmological parameters for each model, and their constraints at la and 2a, obtained from 
an analysis of Union 2.1 SNIa+BAO+CMB+H data sets. 


Parameters 

ACDM 

CPL 

XCPL(I) 

Ao 

N/A 

N/A 

in „ ln -3+2.356xio-3+5.249xio- 3 

+U.DdD X 1U _2.036x10-3-3.643X10- 3 

Ai 

N/A 

N/A 

in fis v 1n -3+2720X10- 3 +3.420X10- 3 
+y.D5 X 1U -6.470x10-3-16.376x10-3 

A2 

N/A 

N/A 

in ,0 v in— 5+2.641 X10 -3 + 5.284X 10 3 
+ /.+J X 1U _ 2 . 644x10-6-3.287X10-6 

CO Q 

- 1.0 

1 09*31 +°- 1118 +°- 1916 

l.UZOl_o 1033 _ 0 .2299 

i noAO+ 0 0453 +°- 1140 

1,uu ™_ 0 . 043 5_0. 094 7 

Wl 

N/A 

, n nO49+ 0 0465 + 0 - 1129 
ru,uu ^-0.0842-0.1.592 

_i_0 nfificc+0.0353+0.1148 

Tu.uuuu_oo638_o.n83 

UJ 2 

N/A 

N/A 

N/A 

&DM, 0 

, n 99fi r+0.0187+0.0475 
TU.^DO_ 0 0241 _ Q 0533 

1 O 99fifi +0.0194+0.0440 
tu '^ oo - 0 .0205-0.0470 

4-0 9971+0.0094+0.0314 

-t-U.ZZ t l-o.oi32-0.0227 

Ho(kms~ 1 Mpc -1 ) 

-1-70 QQf:o + l-5049+4.6282 
' 1 u-yyu ° —1.6108 —4.6937 

- 1-70 9060"’" 0 ' 9572 ' 1 " 2 ' 5967 

-r <u.»uuu_ 1 0491 _ 4 9753 

4-71 n*391+°- 4785 + 2 - 7324 
-1- 1 l.UOZl o.6659-2.1272 

Xtnir, 

567.5023 

559.68457 

557.87849 


Parameters 

XCPL(II) 

DR(1) 

DR(2) 

Ao 

Ai 

A 2 

wo 

Wi 

W 2 

Qdm, 0 

Ho(kms~ 1 Mpc -1 ) 

,n - .n-3+2.356X10-3+5.249X10-3 

-(-U.DcSD X 1U -2.036x10-3-3.643x10-3 

19 2 X ln -3+0. 5X10-3+4.0X10-3 
19.2 X 1U -0.4X10-3-0.5X10-3 

. n .n .n-5+2.641X10-6+5.284X10-6 

-\-Z.43 X 1U -2.644X10- 5 —3.287x10“ 5 

1 nfWl+ 0 - 0453 +0- 1140 

l.UUW_Q Q 435 _Q Q 947 

4-n nfifi^+°-°353+0.1148 

Tu.uuuu_oo638_0.H83 

N/A 

4_0 9971 +0.0094+0.0314 

-i-U.ZZ l l-o. 0132 - 0.0227 

4-71 0*391 +0.4785+2.7324 
' 1 1 - UOZJ -- 0.6659-2.1272 

- 4+3.3243X 10“ 4 +5.4365X 10“ 4 

+0.1 X 1U — 4.7626x10“ 4 — 12.4271 X 10“ 4 

nfi0u in -3 + 107xl0“ 3 +3.5825X10-3 
y.bZ X 1U -1.5925x10-3-4.78x10-3 

+ 1 8 X i q —5 +2.5496X 10“ 5 +4.6683 X 10“ 5 
' -1.7264X10“ 5 -2.6510X10“ 5 

1 0090+0 • 0082+0.0252 
l.UUZU_ 0 0092 _ 0 0346 

+0 ^ X 1 q —5 + 4 -4857X 10 —3 +9.7620 X 10 —5 
+y.O X U _9.7914x10-5-19.801X10“ 5 

fi 0 X i n -9+3.6395xl0 _9 +8.4502xl0 -9 
O.U X 1U -4.4308x10-9-8.9306x10-9 

4-0 9970+0-9085+0.0176 
-f-U.ZZ IU_ 0 0114 _ Q Q269 

4-70 qi +0.7514+3.4634 
- | -' u * yi 1 . 1901 4 .5 339 

- 4+3.3243X 10“ 4 +5.4365X 10“ 4 

+d.I X 1U — 4.7626x10“ 4 — 12.4271x10 — 4 

-3+5.314xl0“3 + ll.4036X10 -3 
+ + 4 X 1U -2.9339x10-3-4.9268x10-3 

+ 1 8 X i q —5 +2.5496X 10“ 5 +4.6683 X 10“ 3 
' ‘ -1.7264x10-5-2.6510x10-5 

1 0090+0 • 0082+0.0252 
l.UUZU_ 0 0092 _ 0 Q346 

+9 ^ X 1 q —5 +4.4857X 10 —5 +9.7620 X 10 —3 
+ y.O X U -9.7914x10-5-19.801x10-5 

fi 0 X i n-9+3.6395Xl0 _9 +8.4502xl0 -9 
O.U X 1U -4.4308x10-9-8.9306x10-9 
4-0 9970+0-9085+0.0176 
-f-u.zz iu_ 0 0114 _ 0 026g 

4-70 qi+0.7514+3.4634 
- | - ,u * yi 1 . 1901 4 .5 339 

Xmin 

558.04200 

556.85491 

556.27088 


TABLE VII. Shows the best fitting parameters today for each model, q 0 , j 0 , I 0 , uj 0 and Udm.o, and their errors at la and 2a, 
obtained from a combination of data. 


Models 

Qo 

3o 

I„ x 10 4 

L0 o 

Qdm, 0 







ACDM 

n rqo9+0.0253+0.0633 
U.oyoz_ 0 . 0331 _0.0700 

-1.00017 +lxl0_O+ i* 10_ ° 

0.0 

-1.0 

0 99CC+0.0187+0.0475 
0241 _ 0 Q533 

CPL 

0 Ol 81 +0.1396+0.2604 
U.Olol _ o. 1454 _ 0 . 326 5 

1 i 47 o+0.2891+0.3985 

1. let i o — 0 . 3092 — 0 . 839 5 

0.0 

1 0981 +0.1118+0.1916 
l.UZCil_ 0 1033 _ Q 22Q9 

4-0 9900+0 0140+0.0440 
+ U.2ZOD_ 0 .0 20 5_ 0 . 0470 

XCPL(I) 

n r Qfi7 +0.0622+0.1622 

U.oyo / -o.0668-0.1362 

1 n8Q9 + 01041 + 0 ' 2047 

i.UO»Z_Q Q 813 _Q 216 g 

r T 141 +23.299+51.964 
° • 11 - 20.093 - 36.098 

i nnzLn+9-9 453 +9- 1140 

x .uueiU — 0 0435 _ 0 og47 

4-0 9971+0.0094+0.0314 
+ U.ZZ <1-0.0132-0.0227 

XCPL(II) 

0 cq«7+0.0621+0.1622 

—0.0669 — 0.1363 

1 n7Q*3+ 0 1047 + 0 ' 2045 

i- u,yo _ 0 . 0835 _ 0 . 2224 

4_r -i i 41 +23.299+51.964 
j. ±etj. _ 20 093 _ 36 09g 

i nrun+ 0 - 0453 +°- 1140 

J-. uuetu _ 0 0435 _ 0 0947 

4-0 9971+0.0094+0.0314 

-fCJ.zz l i_o.oi32-0.0227 

DR(1) 

0 rq 4 c+0.0202+0.0474 

U.jy‘±U_Q Q 25 Q_Q Qg 98 

1 nn8fi+0 • 0262+0.0783 
— 1 .UUOO_q 03 q 6_ o.i203 

_I 9 q 9 nn+3.0693+4.9697 
i-z.yzuu_ 4 3900 _ 12 162 

1 nn9fl+9 • 0082+0.0252 

— l.UUZU_Q QQ 92 _Q Q 346 

4-0 9970+0-0085+0.0176 
+ U.ZZ <u_ 0 0114 - 0 . 0269 

DR(2) 

0 cq 4 r+0.0202+0.0474 

u.uuuu_o.o 2 5o_o,06 9 8 

i nn7fi+9- 0246 +9- 0753 

l.UU/D — 0 . 0299 — 0. 1198 

_ I 9 Q 9 nn+ 3 -0693+4.9697 
TZ.yzuu_ 4 5goo _ 12 162 

1 0090+0-0082+0.0252 
x.uuzu_ 0 0092 _ 0 Q346 

4-0 9970+0-0085+0.0176 
-pU.ZZ U — 0 . 0114 — 0. 0269 


TABLE VIII. Shows the effects of the reconstructions of Iq(z) and w on the parameters Odm, q and j in determined z ranges. 


Model 

XCPL(I) / DR{2) 

XCPL(II)/DR{1 ) 

XCPL(I)/XCPL(II) 

Effect 

Parameter 

Enhancement 
on ^dm 

Suppression 
on Qdm 

Enhancement 
on q 

z 

CO 

I q(I-) 
i q(i+) 

Qdm 

Q 

(+cxd —r -(-3.98) 

(+0.01 -+ +0.1) 
(+0.80 -r +0.85) 

(+oo —r +3.98) 

(-0.01 -> -0.1) 

(+0.80 -> +0.70) 

(+0.1 -!• -0.7) 
(-0.998 -> -1.159) 
(-1.384 -s- -6.24) X 10 -3 
(+1.504 -4 13.976) X 10 -3 
(0.178 -r 0) 

(-0.7 -s- -1.3) 


Model 

DR( 1) 

XCPL(I) 

XCPL(II) 

DRW 

DR( 2) 

XCPL(I) 

Effect 

Parameter 

Enhancement 
on j 

Enhancement 
on j 

Enhancement 
on j 

Suppression 
on j 

Suppression 
on j 

Suppression 
on j 

z 

CO 

Iq X 10+ 3 
Qdm 

j 

(+8 -> 0.295) 
(-1.001 -> -1.002) 
(-76.65 -4 -2.524) 
(+0.591 -> +0.346) 
(-0.9953 -> -1.0) 

( + 1.723 -4 0.045) 
(-0.980 -r -1.002) 
(+17.217 -> +0.968) 
(+0.496 -> +0.258) 
(-0.9968 -> -1.0698) 

(+8 -> 0.045) 
(-0.982 -> -1.002) 
(-153.064 -> -0.320) 
(+0.461 -> +0.241) 
(-0.9821 -4 -1.0698) 

(+0.295 ^ -1.0) 
(-1.0019 -> -1.0021) 
(-2.524 ^ +9.92) 
(+0.349 -> +0.051) 
(-1.0003 ^ -1.0098) 

(+8.0 -> -1.0) 
(-1.0013 -> -1.0021) 
(+19.51 -> -2.09) 
(+0.346 ^ +0.003) 
(-1.0061 -> -1.0098) 

(+8 -> +1.723) 
(-0.9449 -> -0.9619) 
(+77.976 ^ +17.217) 
(+0.496 -> +0.388) 
(-0.9968 -> -1.0176) 


the values of q 0 and j 0 were estimated from a series, in 
where, new variables were defined to avoid the problem of 
divergence, for example, see both Table 2 in [67j and Ta¬ 
bles I, II and III in [681 ] . respectively. These results are 


compatible at ler error with those found in Table IVIII 
Otherwise, constraints on q a and j 0 were established by 
the authors in j69| (see Table 2, series 2D and 3D) from 
a general expression of the BAO modes, and also, em- 
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FIG. 4. (color online) The upper panels display the best reconstructed u(z) along a, their errors at lcr and 2cr for the coupled 
models. These results are consistent with the ACDM model predictions. Similarly, The lower panels show the reconstructed 
evolution of Iq and their errors at 1<t and 2<r in function of 2 for the coupled models and compared with the predictions of the 
uncoupled models. 



log(z) 




iog(4) 



FIG. 5. (color online) The upper panels shows the best reconstructed energy densities at different 2 and their errors at lcr and 
2cr for the coupled models. Here, for 2 > 3.98 and I+, the values of the amplitudes of Qdm are amplified, instead, at this same 
region and for I_ < 0, these amplitudes are suppressed with respect to the uncoupled models. The left below panel shows the 
evolution of Iq as function of 12 dm- For f1 dm > 0.25, a disminution in the concentration of Qdm implies that the amplitude of 
Iq is increased and reduced in the coupled models, instead, for Q.dm < 0.25, it is always increased. Also, the right below panel 
shows the evolution of u as function of SIdm- Here, a disminution of the values of Qdm imply an increasing of the amplitude 
of oj. Moreover, the “o” denotes the best fitting parameters in the present (see Table lYHl) . 


t ing a Taylor expansion and Pade approximations in VI. CONCLUSIONS 

(see Table I, fourth column, fit 3). The results ob¬ 
tained in |69] and 0 at lcr are also consistent with those 
presented in Table IV111 


Now we summarize our main results: 

• An analysis combined of data was performed to break 



































































12 







ACDM 

CPL . ~~ 

Case 1 XCPL 

Case II XCPL 

Case 1 DR . 

Case 2 DR . 

Oo^fOx^PL ^ 


- 0.01 - 0.005 0 0.005 0.01 0.015 


Iq(z) 



'q(z) 


FIG. 6. (color online) The left above panel shows the evolution of q along z. The right panel displays the evolution of u as 
function of Iq. Here, the magnitude of u, |u;| is enhanced from |—0.95| to | — 1.21 when the values of Iq moves from 0 to |0.015|. 
The left and right below panels, depict the evolution of q as function of £Idm and Iq, respectively. Here, the magnitude of q , 
\q\, is amplified from 0.1 to | —1.4| or from 0.4 to | —1.21, when the values of Odm change from 0.6 to 0 or when the values of 
Iq(z) moves from 0 to |0.015|, respectively. Here, the “o” denotes the best fitting parameters in the present (see Table |VTT|) . 
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FIG. 7. (color online) The left above panel shows the evolution of j along z. The right panel displays the evolution of j as 
function of Iq. Here, the magnitude |j|, is amplified from |—0.961 to | — 1.14|, when the values of Iq(z) moves from 0 to 0.005 
and from 0 to —0.001, respectively. The left and right below panels, depict the evolution of j as function of w and SIdm, 
respectively. Here, the magnitude |j|, is increased from |—0.961 to | — 1.121 or from |—0.951 to | —1.3|, when the values of w moves 
from —0.96 to —1.02, and also, when the values of Q.dm change from 0.7 to 0, respectively. The “o” denotes the best fitting 
values of the parameters in the present (see Table Effl . 


the degeneracy among the cosmological parameters of our 
models, allow us to obtain constraints more stringent on 
them. In particular, for the XCPL and DR models, the 
allowed region of its parameters was significantly reduced 
by the inclusion of the CMB data, compared with studies 


of models without the CMB data m[zi. This implies 
that higher redshift may be able to discriminate between 
these models. 

• In the DR model, a novel reconstruction for to was 
proposed whose best fitted value is closed to —1, and 
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has the property of avoiding divergences in a distant fu¬ 
ture z —1. This result is consistent with the value 
predicted by the ACDM model at lcr error. Lisewise, 
within this scenario, a finite value for oj has been ob¬ 
tained from the past to the future, mainly, the asymp¬ 
totic values are: u>(z) = LO 2 Z 2 for z 1, ui(z) ~ u> 0 for 
z«l and uj(z) ~ uj 0 for z —> —1. Therefore, a better 
physical description of the dynamical evolution of DE is 
performed by the DR model, which should be used to 
explore the properties of DE. 

• Currently, a phase of accelerated expansion is the situa¬ 
tion revealed by all our models about the universe. The 
big-rip problem is not forecasted by the DR model, and 
hence, this scenario should be considered to study the 
ultimate destinity of the universe. Likewise, from the 
coupled models found that the values of the amplitudes 
of the parameter q are not significantly affected neither 
by values of Iq nor by the values of w (see left above 
panel in Figure ©. 

• The values of the amplitudes of £Ide (see upper panels 
in Figure [5]) are not significantly modified by the recons¬ 
tructions of Iq and w, respectively, nevertheless, they are 
definitely positive. This requirement implies that ui must 
be always negative in all the cosmic stages of the universe 
(see upper panels in Figure [I]). 

• The coupled models are strongly favored by the obser¬ 
vational data having a preference for jo < —1, and hence, 
they represent a slight deviation from the value predicted 


by the ACDM model (see Table IVIII) . 

• The behaviours qualitatively presented here show that 
the graph of j has more possibility in discriminating the 
different coupled DE models, and therefore, j could be 
used to distinguish them (see left upper panel in Figure 

ED. 

• The physical effects generated by the magnitudes of 
Iq and ui on the cosmological parameters could be un¬ 
derstood, so: An energy transfer from DE to DM or 
vice versa inserts energy into one of the fluids, and deter¬ 
mines an increase of the energy density on one of them, 
which increases the Hubble parameter inducing a slight 
expansion of the universe, and to recover equilibrium of 
the system Iq and ui leads to an enhancement or sup¬ 
pression on the amplitudes and shapes of SI dm, q and j 
with respect to uncoupled models, in determined redshift 
ranges. 

In a forthcoming paper we will extend our study by ap¬ 
plying cosmological perturbation theory on the coupled 
models, using data of linear matter power spectrum, 
weak lensing potential, integrated Sach-Wolfe, growth 
rate, and other. They will allow us to calculate of how 
the magnitudes of Iq and u ; operate on the amplitudes 
of S 1dm, q and j, respectively. From which, we may con¬ 
clude if our DR model can emerge as an alternative to 
the ACDM model. This will be the purpose of our future 
work. 

Appendix: Integrals I n (z) and I n (x ) 
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